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Abstract 
A  nonlinear  boundary  value  problem  is  considered  for 
the  axisyininetric  buckling  of  thin  spherical  shells  subjected 
to  uniform  external  pressure.  The  uniformly  compressed 
spherical  state  is  a  solution  of  this  problem  for  all  values 
of  the  pressure.  Me   prove,  using  Poincare's  method,  that 
for  pressures  sufficiently  near  each  simple  eigenvalue  of 
the  linearized  shell  buckling  theory,  there  is  another 
(buckled)  solution  of  the  nonlinear  problem.   A  convergent 
perturbation  expansion  is  used  to  analyze  the  buckled  solutions 
near  the  eigenvalues.   For  a  limited  range  of  caps,  we  also 
prove  that  one  or  three  buckled  solutions  bifurcate  from  the 
multiple  (double)  eigenvalues  depending  on  their  order.   The 
existence  of  a  "lovfest"  intermediate  buckling  is  established 
and  precise  upper  and  lov/er  bounds  are  given  on  its  magnitude. 

1.  Introduction 

The  surface  of  a  thin  elastic  spherical  cap  is  subjected 
to  a  uniform  pressure,  p,  which  is  directed  towards  the  cap's 
center  of  curvature.   It  is  a  vrell  known  experimental  result, 
see  e.g,  [l] ,  that  as  the  pressure  increases  from  zero  the 
cap  deforms  only  slightly  from  the  spherical  shape  until  a 
critical  pressure  p  =  p.  is  reached.   Then  the  cap  suddenly 
jumps,  with  relatively  large  deflections,  into  a  non-spherical 
shape  which  we  call  the  buckled  state.   The  fundamental  problem 
of  shell  buckling  is  to  determine  the  mechanism  v/aich  initiates 
the  jumping  and  to  obtain  estimates  of  p. . 
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Previous  investigators  I  have  assumed,  as  an  approxima- 
tion to  experimental  conditions,  that  the  cap's  edge  is  rigidly 
damped,  i,e,  the  displacement  and  change  of  slope  are  zero. 
In  this  paper  edge  conditions  are  considered  for  which  the 
spherical  shape  (the  unbuckled  solution)  is  a  possible  solution 
of  the  nonlinear  problem  for  all  pressures.  V/e  refer  to  these 
as  bifurcation  buckling  problems  since  other  (buckled)  solu- 
tions of  the  nonlinear  problem  may  branch  from  the  unbuckled 
solution."**  For  example,  a  bifurcation  buckling  problem,  which 
we  call  Problem  B,  is  obtained  if  the  following  conditions  are 
specified  on  the  edge  of  the  cap:   no  rotation  (clamped); 
zero  transverse  shear  force,  i.Cc  the  edge  is  free  to  move 
normal  to  the  spherical  surface;  and  the  meridional  membrane 
stress  is  prescribed  so  that  it  is  in  equilibrium  with  the 
applied  surface  pressure.   Other  bifurcation  problems  are 
obtaJ.ned,  for  example,  by  replacing  the  condition  on  the 
meridional  stress  in  Problem  B  with  a  corresponding  one  on  the 
meridional  membrane  displacement,  or  by  permitting  the  cap  to 
freely  rotate  instead  of  clamping  it.   The  bifurcation  problems 
are  precisely  formulated  in  Section  2,   In  Sections  3,  I4.  and  5 
only  Problem  B  is  considered.   However  analogous  results  can 
be  established  for  other  bifurcation  problems.   Some  of  these 
are  contained  in  the  final  section  of  the  paper.  We  consider 
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1  See  e,g,  [l-7j  and  references  contained  therein. 

'"■  For  the  rigidly  clamped  cap,  which  v/e  refer  to  as  relaxation 
buckling,  the  spherical  shape  is  a  solution  if  and  only  if 
p  =  0.  ... 
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only  axisyranetric  deformations  of  the  cap. 

Bifurcation  buckling  problems  for  spherical  caps  were 
first  considered  in  an  approximate  form  in  [8] ,   An  equivalent 
of  Problem  B  was  treated  in  [5]  where  the  linearized  buckling 
theory  was  partially  analyzed  and  approximate  solutions  of 
the  nonlinear  problem  were  obtained.   The  linear  buckling 
theory  had  been  previously  discussed.'"' 

The  precise  knovjledge  of  the  unbuckled  state  for  the 
bifurcation  problems,  permits  us  to  rigorously  establish 
certain  properties  of  the  solution.   In  Section  3  we  prove 
that  for  all  P  sufficiently  near  each  simple  eigenvalue  of 
the  linearized  shell  buckling  theory,  a  solution  of  the  non- 
linear problem  exists.   Only  one  solution  of  the  nonlinear 
problem  branches  from  each  simple  eigenvalue.   A  related  result 
is  established  in  Section  I4.  for  the  other  eigenvalues,  which 
are  all  double,  and  for  a  limited  range  of  caps.   However,  we 
find  the  surprising  result  that  one  or  three  solutions  of  the 
nonlinear  problem  bifurcate  from  the  double  eigenvalues 
depending  on  their  order.   The  conjectured  load-deformation 
curves  for  simple  eigenvalues  are  sketched  in  Fig.  1.   A 
perturbation  expansion,  which  is  valid  near  each  eigenvalue, 
is  used  to  prove  that  the  curves  in  Fig.  1  have  the  form  of 
the  solid  portion. 

Our  method  of  analysis  is  the  Poincare  bifurcation 
theory  used  to  prove  the  existence  of  periodic  solutions  of 


Privately  comrr.unicated  to  the  author  by  \-l ,    Squires,  1958. 

Here  P  is  a  dimensionless  parameter  proportional  to  p,  see 
Eq.  (2.2a)  below. 
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initial  value  problems.   Previous  applications  of  this  tech- 
nique to  boundary  value  problems  are  given  in  [9  ,  lO]  v;here 
related  buckling  problems  for  columns  and  circular  plates  are 
studied.   In  Section  3»  where  the  simple  eigenvalues  are 
studied,  the  procedures  employed  are  closely  related  to  those 
given  in  [lO]  .   Modifications  are  made  in  Section  l|.  to  inves- 
tigate the  double  eigenvalue  case. 

Friedrichs  [ll]  proposed  an  energy  mechanism  and  intro- 
duced the  concept  of  an  intermediate  buckling  load  to  explain 
the  experimentally  observed  juraping  of  complete  spheres  from 
an  unbuckled  to  a  buckled  state.   Modifications  and  extensions 
of  these  ideas  were  subsequently  proposed  in  [l2]  .   As  applied 
to  the  bifurcation  buckling  of  caps,  see  e.g.  Pig.  1,  sev- 
eral intermediate  buckling  loads  pA^   corresponding  to  differ- 
ent branches  of  the  solution  may  exist.   For  a  given  branch, 
say  bifurcating  from  the  eigenvalue  P  ,  Pj:    is  defined  as 

a  load  in  the  interval  P;"^  <  P..   <P   such  that  for  all  P  in 
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Pt   _  P  "^  ^M   ^'^®  unbuckled  state  has  less  potential  energy 
than  the  corresponding  buckled  state  and  conversely  for  P  in 
pA   <  P  <  P  .   In  Section  5  v;e  establish,  assuming  that  the 
potential  energy  has  a  minimum  for  every  finite  P  and  P  >  0, 
the  existence  of  a  lov/est  intermediate  buckling  load  P.,,  see 
Fig.  1.   It  is  shown  to  be  bounded  from  belol^r  by  the  lowest 
buckling  load  of  an  "equivalent"  flat  circular  plate  buckling 
problem  and  bounded  from  above  by  the  lowest  eigenvalue,  ?,    of 
the  linearized  shell  buckling  theory.   More  accurate  upper 
bounds,  which  are  considerably  less  than  _P,  are  rigorously 
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obtained  in  Section  5  by  a  nininization  procedure.   These 
bounds  bracket  P.,  for  a  limited  range  of  caps.   It  seems 
likely  that  v>rith  suitable  modifications  corresponding  esti- 
mates of  intermediate  buckling  loads  can  be  obtained  for  the 
unsymmetric  buckling  of  spherical  caps  and  other  shell  bifur- 
cation buckling  problems.   Some  of  the  results  of  Section  5 
are  closely  related  to  ones  previously  announced  by 
Vorovich  [iSj  • 

2.  Formulation  of  the  Boundary  Value  Problem 

The  elastic  spherical  cap  is  of  thickness  2h  and  radius 
R  and  has  a  small  angle  of  opening  2/\  ;  see  Fig.  2  for  the 
shell  geometry,  V/e  consider  the  axisymmetric  deformations  of 
the  cap  that  result  from  a  uniform  and  inwardly  directed  pres- 
sure p.   The  non-vanishing  mid-surface  displacements,  u  and  w 
which  are  in  the  meridional  and  normal  directions  to  the  shells 
mid-surface  are  therefore  functions  only  of  the  polar  angle  9, 
Both  w  and  p  are  counted  positive  when  directed  towards  the 
center  of  curvature.   Nonlinear  differential  equations  which 
describe  the  small  finite  axisymmetric  deformations  of  thin 
spherical  caps  have  been  derived  by  several  authors,  e.g. 
[l,2,ll|.-l6]  .   These  equations  may  be  written  as, 

(2.1a)     Gf(x)  +Xf(x)  =  pg(x)[f(x)  +1]  , 

(2.1b)      Gg(x)  =  -P[f^(x)  +2f(x)J  .  :. 
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Here  we  have  used  tlie  dimensionless  variables: 

X  =  o//v  ,   f(x)  =  ^  ^'^^   ,   p  =  (A^/C)(R/h)  , 


(2.2a) 


P  =  (R/h)^(p/2EC)  ,   X  =  PP  ,   C^  =  2/3(1  -v^)  , 


and  the  linear  differential  operator  G  which  is  defined  by 

G(})(x)  =  x-^[x^^'(x)]'  , 
where  a  prime  indicates  differentiation  with  respect  to  x. 
Here  E  is  Young's  modulus  and  v  is  Poisson's  ratio.'"* 

The  "excess"  stress  function  g(x)  in  (2.1)  is  defined 
such  that  the  meridional  and  circumferential  membrane  stresses 
Oq  and  oi  and  the  corresponding  dimensionless  stresses  /^fi  and 
r,  are  given  by 

ZqM    ^    (R/h)(2/EC)0Q(9)  =  g(x)  -P  , 
(2.2b) 

^(x)  =  (R/h)(2/EC)o^(9)  =  [xg(x)  -xP]'  . 

The  outer  surface  bending  stresses  o^,  and  oi  ,  and  the  corre- 
sponding dimensionless  stresses  ]r.  and  ^i  are  given  in  terms 
of  f(x)  by 

ZlM    =  (R/h)(2/3S;c2)o°(9)  =  xf'(x)  +  (1+ v)f(x)  , 
(2.2c) 

^(x)  H  (R/h)(2/3Ec2)o^(9)  =  vxf'(x)  +  (1+  v)f(x)  . 

The  independent  variable  f(x)  is  related  to  the  slope  of 
the  deformed  middle  surface  of  the  cap  with  respect  to  the 
initial  spherical  shape.   Thus  if  f(x)  =  0  for  a  given  defor- 
mation the  deformed  middle  surface  is  also  spherical.   We 

The  independent  variables  a(x)  and  y(x)    employed  in  previous 
papers  [2,3]    are  related  to  the  present  variables  by, 

a(x)  =x[f(x)+lj  ,   y(x)  =x£g(x)-Pj  . 
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refer  to  P,    which  is  defined  in  (2.2a),  as  the  georaetrical 
pararaeter  and  P  and  \   as  either  pressure  or  loading  parameters, 

To  complete  the  formulation  conditions  at  the  center 
X  =  0  and  the  ed^e  x  =  1  are  required.   From  the  symmetry  of 
the  deformation  and  the  regularity  of  the  membrane  and  bend- 
ing stresses  at  the  origin  we  obtain  with  aid  of  (2ol)  that, 

(2.3)  f'(0)  =  g'(0)  =  0  . 

The  edge  of  the  cap  is  restrained  from  rotating  so  that, 

(2.l|a)  f(l)  =  0  . 

In  addition,  ^^fe  assume  that  at  x  =  1  the  transverse  shear 
force  vanishes  and  ^^(1)  is  specified  so  that  the  cap  is  in 
equilibrium  v/ith  the  applied  pressure.   This  yields  the  bound- 
ary condition 

(2.i;b)  g(l)  =  0  . 

The  bifurcation  buckling  problem.  Problem  B,  is  defined 
as  the  boundary  value  problem  consisting  of  the  differential 
equations  (2.1)  and  the  boundary  conditions  (2.3)  and  {2,[\.)» 

Other  boundary  conditions  can  be  specified  at  x  =  1  -  ■ 
to  yield  bifurcation  buckling  problems  for  clamped  caps,  e.g. 

(2.5)     f(l)  =  0  ,   g'(l)  +  (1-  v)s(l)  =  0  .  , 
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These  conditions  imply  that  on  the  edge  the  transverse  shear 
force  vanishes  and  the  meridional  displacement  u  (or  equiva- 
lently  the  horizontal  displacement)  is  proportional  to  P   The 
proportionality  constant  is  determined  so  that  the  spherical 
form  is  a  possible  solution  for  all  P. 

3,  The  Existence  of  Buckled  Solutions. 

A  solution  of  Problem  B  that  is  valid  for  all  finite  P 
and  r  is  the  unbuckled  solution,  f  (x)  ?  g{y-)    =  0.   This  corre- 
sponds to  a  state  of  uniform  compression  in  v;hich 


^  (x)  =  ^^  (x)  =  -P  , 


and  the  deformed  middle  surface  remains  spherical. 

The  existence  of  buckled  states  will  now  be  established 
using  Poincare's  method.   Specifically,  we  prove  in  this  sec- 
tion that  for  every  -•'  >  0  and  for  every  positive  integer  n 
there  exists  a  buckled  solution  when  P  is  in  a  sufficiently 
small  and  full  interval  about  the  n-th  simple  eigenvalue, 
P  ,  of  the  linearized  shell  buckling  theory.   This  is  in  con- 
trast to  the  buckled  circular  plate  [lOj  where  the  load  must 
be  slightly  preater  than  each  eigenvalue  of  the  corresponding 
linear  buckling  theory.   I7e  also  show  that  for  P  sufficiently 

near  P   the  solution  has  n-1  internal  nodes.   In  the  follow- 
n 

ing  section  the  solutions  of  Problem  B  near  the  other  eigen- 
values of  the  linearized  theory,  which  are  all  double,  are 
studied. 


'")  .f.  "■,■)! 


"A   ':ic;;    ;:.lJ.bv    ox    C';.";j    :'   iixo :..:•;'.■  !       :'   ij;/io' 


;..i-i     ,:;  -   ;-):^    ^    {< 


n.r 


I  ■ 


-oar    ^Id^    a.j  ■^veaa    ©-•!   , '^I.C^-oi:  i  1  n':^'.;:::       ■;.'-•-■'•      •   ;:':-^;-'   ^:^^   ?:^''M^-f^ 

:;    '■.::•  J;  y  J*  fj  I  avJciaoq    ■^'i6V&    l.^"    b;':::    0    f^  -     :-;;::-e    ■•X'-.i.    V  si.-n    n.,NJ;;t 

\  [.y       '  .Oils,    ■'         ■                <                               '                               ■ 

-woIX  "      ,',->•■'-■,-    Camacfax  I-a  ear*  r-'-. -.,        ■-,„   fd;f   ^'i 


To  employ  Poincare's  method  parameters  e  and  "'and  new 
independent  variables  y(x)  and  z(x)  are  defined  by, 


(3.1a)    e  =   lim   [xf(x)]'   ,    6  =  e"""-  lin   [xg(x)]  '  , 
X  — ^  0  X  — ^  0 


(3.1b)    y(x)  =  e"-'-xf(x)     ,    z(x)  e  e"-^xg(x)  . 


The  differential  equations  (2.1)  of  Problem  B  are  then  given 
by, 

(3.2a)   Hy(x)  +Ppy(x)  =  oz  (x)  [l  +  ey  (x)/x]  , 
(3.2b)   Hz(x)  =  -p[ey2(x)/x +2y(x)]  , 

where  H  is  the  linear  differential  operator, 


(3.2c)         hMx)  =  jx-^[xMx)]'r 


The  initial  value  problem,  ^  ,  is  defined  by  the 
differential  equations  (3«2)  and  the  initial  conditions: 

(3.3a)        y(0)  =  0   ,    y' (0)  =  1  , 

(3.3b)        z(0)  =  0   ,    z' (0)  =  b. 

For  all  finite  P,s,  "^j  and  P  >  0,  Problem  *'.'■    has  a  unique 
solution  in  the  interval  0  _<  x  _<  1  which  is  analytic  in 


■•>if    .•■  i'/j 


}«.^;    3 


'V"?  "        .'if.:  ''-I-  • 


'' '  i  T  c  'j . :  ■■  o  \   ■<■'  d  i."  r-  ff:  ^ 


!!>,-.'    :^^r'ii ''Ifsv    'v-fn    ■:■'):.    ;.in;i    i')-,.    b  alJa.i '"i.^v    d'.iotr.o'ioi-ifl' 


\  X  I  .'^ 


'''^  ^   i.^:.)v:       "I'l*';) 


';i-^V.i:-' 


.t.O  I.  y  ...     '   --     ■  J..  ±J-:-:>       ^t  :  i. 


.&  !(:■):. 


(  '^)  -■■         -r    i  -/.'ii':"! 


<  •;.     ■,  f 


^^^    :m- 


\x).jH 


*^r;:i'    <' 


;     ■iV'.r    e.r    ii.    ^?'l«;-fw 


M(^o^::]^""x 


^^■.s*:) 


.    I   =    (0)'\ 


'J 


t'G)'i  (i-^C,  »C.) 


0   =    (0)^ 


(ciC.r 


.•0 
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P,  e  and  6,      This  result  can  be  proved  by  tae  sarae  nethods 
used  in  [lO]  to  establish  an  analogous  result  for  the  buckling 
of  circular  plates  and  hence  the  proof  is  not  given  here.   For 
fixed  p  >  0  values  of  the  parameters  P,  e  and  rf  are  sought 
such  that  the  solution,  [y(x;P,  e, 6)  >  z(xjP,e,o)J  off-*''  satisfies 
the  boundary  conditions 

(3.L^)       y(l;P,e,o)  =  z(l;P,e,.5)  =  0  . 

If  such  a  choice  of  parameters  is  possible  it  then  follows 
from  (3.1b)  that  f(x)  =  ^y(^^P>g'-^)  and  g(x)  =  ez(x;P.e.6) 

is  a  solution  of  B. 

A  solution  of  >-^- vjhich  satisfies  the  boundary  conditions 
(3,14.)  is  obtained  by  choosing  the  special  parameter  value, 
e  =  0  in  (3.2) .   The  resulting  initial  value  problem  for  the 
functions  y  =  y (x-,P,0,  •^)  ,  z  =  z(x;P,0,:-')  has  solutions  which 
satisfy  (3.I4-)  if  and  only  if 


(3.5a)   P(.)  =  P^(p)  =  2r/^  +  .g/p  >  2  /2 


2  .  2 

'P   ^   c  yd.  j 

n  =  1,2,..., 


(3.5b)   6(r)  =  ^^ip)   -  20/^1  J 


where  '-  is  the  n-th  zero  of  the  Beasel  function  J,  (x).   The 
quantities  P„(f")  are  the  eigenvalues  (or  buckling  loads)  of 
the  linearized  buckling  theory  and  they  are  simple  eigenvalues 

^^      P^     Pmr.  a^<i  P  ^     ^J^rnJ      =  ^rr^^f^J      WherO, 


:; 


J      -■       \  :•%-'■:        ^  -i    ,    •■        ■  ■  t.  -     ■      '■  r      -   '    V^  V   '^»  <-.  ^ 

.  3  ■'■  u   fjoi ;'j;Iao    n   ax. 


<     O 


"f(  ft       '   '  1  i 


10   (ebsoX  snxX^IojJd  io)   eeiJXBvnss-t©  s"^^  sis   (■'-)^'i  a3i.;JxJnflyp 

61.0  ,        .  ,^^.  ^    ^^j_^ 
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(3.6)    P^^   ^  2  ^  <c'  c>  ,   m,n  =  1,2,...,    m  7^  n  . 


The  solutions  are  the  eigenfunctions  (or  buckling  modes)  of 
the  linearized  theory  and  for  simple  eigenvalues  are  given  by""' 


(3.7) 


'  n  —  x^i^,»w*. 


v;here 


z(x;P^,0,cS^)  =  z^"^x)  =  A^J^C-^x)  , 


A  =  2'A^j'(0)  . 
n     '^  n  1 


The  corresponding  solutions  of  Problem  B  reduce  to  the  unbuck- 
led solution. 

Thus  for  each  finite  p  >  0  there  are  a  denumerably 
infinite  number  of  roots  1?^'^*^^*  ^  ~  1»2,...  of  O.lj.). 
According  to  the  implicit  function  there  are  other  roots  of 
(3.I4.)  near  i^^f^f^'.^^    if  the  appropriate  Jacobian 


(3.6)   J  =  z["^l)y(^)(l)-z(^^l)y(")(l)  ^0  ,     n  =  1,2,...  . 


h 


Here  we  have  used  the  notation 


For  simplicity,  the  explicit  dependence  of  the  solutions  on 
p  is  suppressed  here  and  in  the  remainder  of  the  section. 


i.l 


a  \  ?t 


sxr-. 


p  *  V      t 


rtni 


f   -■»*•■*  -.V  ^.  . 


(v.c) 


m 


in'f  •'l'-i^! 


n/j-'^iocol.   scJ-fiiTiro'iqQB   8.do    'ii    ;!  _^,'.  i..  ■':.t^^;-l]    ijjorj    (.j-Lt) 


_,S,I:=ri  ,    0^    {l)^'':xU.)'-''U-{n'''U{i)^''\^   ^l      (3.0 


noi;ii?Jon  ©x"J^  t'j'ajj  avari   ow  gtoH 


no   c 
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(3.9)    T-   y(x;P,e,.-',) 


d'b 


=   j]^'  (x)  ,  etc 
o 


n=  1,2,...   . 


F=P 


n 


6=0 


n 


Since  y  and  z   are  analytic  in  P  and  <3,  the  initial  value 

problems  (variational  problems)  satisfied  by  y- ,  y, ,  Zp 

and  z^  can  be  obtained  by  formal  differentiation  of  -^  with 
b 

respect  to  P  and  (5.      Thus  we  obtain  from  (3.2)  and  (3.3) 
using  (3-7)  and  (3-9), 


(3.10a)  Hy(^^P^py(^)=p[z(^).y(^)]  ,  y^^^O)  =  y^^'^^'(0)=o) 


(3.10b)  Hz^=^^  =  -2^7^''^  , 


(3.11a)  Hy^^^P^ry^)  =  P3^)  , 


z(^)(C)  =4^^0)=0;! 


yl'^No)  =y^T'''(0)=0, 


Sn=l,2,.,. 


( 


(3aib)  Hzi""^  = 


I? 


(_")(0)  =0,zHo)=l], 


The  explicit  solution  of  (3cll)  (not  presented)  yields, 
at  X  =  1,  the  relations 


'V  0  if  T'  ^  r. 


(3.12)   zi"^l)  =2.o-ly(-"^l)  s 
5        no' 


mn 


f  in, n=l,2,...,  TTi/^n • 


=  0  if  p  =  P 


mn 


'  .;ic'  .0'-,  - 


in) 


\^    - 


f- 


■J  J.vi    ■■    ■!•>   rox  J, :;,;:.;!■.  A' ;-  ;.'':  ;  ■!> 


-■i.i  iv;3'j 


I '^i   :;   J..  ..»■    \;   ^^)alf\ 


;,o-:o)- 


■  ■,  ^  ^  -  -  . 


■:-v;H    ■  .G  J.f . .-  ) 


.o-=(.)'^^?..^   r.-r^'-U 


..n^v*^:- 


:i«fo,'^^^3,o-  io)^^^;', 


Vi-:; 


.',t-. 


v-w. 


ii  u--:..j 


i;-"  -  '-'Kr  (j:.c.5.) 


n.-^  \^1s    Ok 


.a\m  ,  .  , .  iS,I=n<,oi  , 


m. 


Ix   0  = 
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Therefore    it   follov/s   from    (3.8)    and    (3.12)    that 


(3.13)  J  =  V^''Ul)\^-\i''Ul)   -z^^^l) 


n  -^P 


Hence,  J  =  0  if  P  =  T   and  P  =  P   and  nothing  can  be  concluded 
concerning  tne  existence  of  neighboring  roots  to 
IP  (P   )*0,-'^  (r'  )1  .   This  case  is  considered  in  the  follov/ing 
section.   To  prove  that  J  7^  0  if  p  7^  I'^w.r,'  P  7^  ?„  ^e  assume 
the  contrary,  i.e^  let. 


'  '  ■  mn'  '    '       n 


(3. Ill)        y^^'^l)  =  2-\z^^^l)  =  a 

where  a  is  an  arbitrary  real  constant.   Then,[yp   (x),Zp   (x)^ 
is  a  solution  of  (3.10)  that  satisfies  the  boundary  conditions 
(3»1^).   A  solution  of  this  boundary  value  problem  exists 
if  and  only  if  the  in homogeneous  terms  are  orthogonal  to  every 
solution  of  the  homogeneous  adjoint  problem.   This  leads  to 
the  condition  that  for  all  a 

,1 

/   J^(6'^x)xdx  =  0  , 
"0 

which  is  impossible  and  hence  J  7^  0. 

Hence,  by  the  implicit  function  theorem  and  the  analyti- 
city  of  y  and  z  in  the  parameters,  (3.^4-)  can  be  uniquely 
solved  for  P  and  6  as  analytic  functions  of  e  in  some  suf- 
ficiently small  neighborhood  of  each  root  {?    y  0 ,r:   ]  ,   p  ^  p 
The  solutions  are  the  analytic  functions 


iiii-];}    (Si-J:)    '-rif^    (O.l)    l'O-^t    ?  joI 


r,->'i     -i 


s'lo'ioT-.TiT 


_    '  r 


.  '■'  .1 


n^^-x- 


'  1.  V-,;.- 


oiriua^ii   ^w  ^^^s   i-   'i    (.,.'   -^  ~; 


li 


\- 


l-    ^,'3,  ;;    ovo'i',:; 


■''''"''  '';T';    ■'•"i"'^^ 
,  y;  C'  i  :i ; j  ?  (i 


':)••' 


-^ 


■•■)^'';.7 


(ii-'-^) 


,,.  r  .^ ,  •  r  ^  ■ 


.    0   -   xi:vc(z^'OfV 


A 


.0  ^;  l-   eorio.i   bnsi   a.L'dicsoqrfii    si  jiuldw 
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(3.15)      P   =  P(n)(^)»       '^=    '^(n)^^^'       '^1    "^    ^n'      "  "   1,2,...    , 


v;hich   satisfy  the    conditions 


P(n)(°)    =   Pn    '      ^n)(^)    =  ^n 


For  sufficiently  small  e 


f(x;e)  =  ex~-'-y^(xie)  =  £x"-'-y(x;P  ^^^  (  e)  ,  £,.5(^)  (  e))  , 
(3.16) 

.-1_   .„..N -1 


g 


(x;e)  =  ex"  z^{x',t)    =  ex"  z  (x;P  ^^^  (  e)  ,  e,  ^^^^  (  e)) 


are  solutions  of  B. 

The  solutions  of  B  near  e  =  0  are  now  considered.   Since 

P/  ^(c)»  '^/  \(e)>  Y  (>^;e)  and  z  (x:e)  are  analytic  functions 
(n)  ^  '  '  -  (n)  ^  '  '  -'n   '        n  "^ 

there,  they  have  convergent  expansions  in  some  interval 
I e|  <  e°  given  by: 


n 


CD    /„x   .  00 


(3.17) 

y^U;^)    =  y^"^x)+T^y[^^x)£\  z^(x;e)  =2^"^x)-f7~-z["^(x)£^ 

The  expansion  coefficients  are  determined  in  the  usual  way  by 
substituting  (3«17)  into  (3«2)  and  {3>»U.)    ^^^   equating  the 
coefficients  of  each  power  of  e.   Each  of  the  resulting 
system  of  linear  boundary  value  problems  for  the  coefficients 

[yi   ,  z^  ']  has  a  solution  if  and  only  if  the  appropriate 


-:  ) 


•\       :'.  r'  -r      ^, 


V  i^i: HI.   frj  Li-r 


^1  \ 


r";., ,:     ^'i,'    :i--J.  0  -  :    '^iJH     '..:}'■■ 


r.-.-'^V 


)..-.;  .^.(o.^'iv--^-' 


',      e  -• '  /    ■ 
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■  I 
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"  l-;x),^,: 


\i.\  v£^w   Ig; 
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orth.ogonality  condition  is  satisfied.   For  the  first  two  of 
these  problems  with  i  =  1,2  the  orthogonality  conditions 
yield 

1 

1 

(3.18b)   P^'^)  =  ^^-^  ;'   [2y|^)(x)+^-lzj^)(x)Jj?(::x)dx 

^  n'^S^  n^^   0 


(3.18a)   P,^"^  =   ,  .  ^\ /   J^('-'x)dx  >  0  , 


2j;(0)o^p(n)   1 


i^-n)    ^ 


(n) 


y|'^Mx)J^('i^^x)xdx  . 


Thus  it  follows  from  the  first  of  (3.17)  and  (3.l8a)  that  for 
each  n  solutions  of  B  must  exist  for  all  P  in  some  small 
interval  about  P^. 

For  sufficiently  small  e  v-xe  have  from  (3.17) 

y^(x;e)  =  y^'^^x)  +0(e)  . 

Hence  in  the  interval  0  <  x  <  1,  y   has  the  same  number  of 

(n)  U  T       •      T 

zeros  as  y    ,  or   y    has  n-1  simple  zeros. 

Approximations  of  the  lov/er  buckling  load,  P/  w, 

for  each  branch  may  be  obtained  by  truncating  the  series  in 

the  first  of  (3.17)  and  from  the  condition  that 

^^n)(^) 

*-T-t =  0  at  the  lower  buckling  load.   These  approximations 

may  be  valid  only  if  the  expansions  (3.17)  converge  for  suf- 
ficiently lar^je  e.   For  example  if  2  terms  are  retained  in  the 
series  so  that 


f  \i 


■ic   qwj   3r'jdl   er::.t  •T-y:      ,-...s^^z]:iH-    bj    r:  [,t.:i3rioo  Y^-cIi^nori,.,  ^rr, 


snorrf-Jl.nc.;    7:t  i:  :ivnv.3o:'c+'5:c    o/;j    r^,. 


L-:.:s'-   BrrreXdoiq    oi'&rj? 


,.    0   <    -'  ■-  ■  ^'^■-  ^  'i'- 


(3or^n 


I. 


f  ,-''>"' 


■t.t) 


^^. 


}:^.(x}^^^       \ 


(T.'«'~}    .^ao-:^^    n^-tiil    V5u<    >    X  [.iu,.:^   ',  '^■iT'ii^  ■  oi."  ^uf-    to'i 
lo  'ioG:uJi'  '.-ni^a   oii.;   e;.,.ti  ^^    ,.[  .->   x  >  0  Imv^: 'cfai    aricf  fu   eoctsH 

v-tBxi*  norctibnoo   sd^  xno'il  bnh    \VX.£)   lo  o  .  . 

"'ma   lol    oa'"i'-iViiv'0    (Vr?'")    Rao  j:-;n.G^  "■'•■•    nri  f   ■>  ^    --f.---^    r-,  r  rnT.'    -,,;    ■-■-. 
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P  (e)  ^?    +P,^'^^e  +P^'^^e^ 
n      n   1      2 


^•^^^  p(n)2 

(n)L   n  Kp(n) 


I4.,  The  Existence  of  Buckled  Solutions  Near  Double  Eigenvalues 

We  now  consider  P=  ■'^   and  the  double  eigenvalues, 

'     mn 

(li.l)  P  =  P  {'      )  =  P  (-'^   )  =.J^+-^  ,   m,n  =  1,2,...,  n  ?^  n  . 
^^   '      m  inn     n  mn   ~  m   n     '       '       ' 

The  solutions  of  the  initial  value  problem  (3.2),  (3.3)  with 

e  =  0,  p  =   r-        and  P  =  P  (/-'  ),  that  satisfy  the  boundary  con- 
*  '    '  mn  m  •  mn  "^  •' 

ditions  (3.i|)  are  the  eigenfunct ions'" 

y(^'Pm(^r.m)'°'^)=y^'"^(-'''')  =  (-mn/-m^)Jl('V)^(-nm/WJl^ 
{i4..2a) 

z(x;P  (p   ),0,r^)  =  z^™^^(x;c)  sa  J,('-^x)-^a   J,(^^->x), 
ra  *  mn  '    *    '  ""  *'^"'   mn  1  m     nm  1^  n   ' 

where    cb  is   an  arbitrary   real   number   and 

e     ^2i/2,.A>    , 

mn  n'    ra 

(14-. 2b)  ,  /p 

^  J1(0)(-'2.,^2)    >-  nm  J     . 

Here,  and  in  all  subsequent  equations  in  this  section, 
the  indicies  m,n  =  1,2,...  ,  m  ^   n,  and  we  shall  hereafter 
omit  explicit  reference  to  this. 


L        .    -I  .      .     ■         ;i 
in), J 


r  ~ 


acrid" 


<•.'  ...■ 


*    (nU,,     a:"  J(k; 


.  2e;."Isvn^;;  io   v».rdjjc^:    ■^if^-i   bciB       '''■    -  -■    -.oMiinTOo   von   v^w 

t  .  ■  .1,  i  i*i  I.  ->'J  .'.  4.-1*1  4.1  ' 

n.'-^  ■  .III  nr.i  ' 


*       ..  ■ 
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The  double  eigenvalue  case  may  be  of  especial  physical 
significance.   Consider  the  variation  with  p  of  _P(.'^)  v;hich  is 
defined  by 


(U.3)  P(P)  =  min  P^(p) 

n 


There  are  local  minimums  of  P([^)  at  P  =  2'''     ,  p   =  2 


=  ?3/2  ,„  ^  pl/2^,2 

n 
dP 


n  =  1,2,...  .   At  p  =  P   .T  -r-z,   is  discontinuous   or  there 

'    ' n,n+l  dp 

are  "peaks"  in  the  curve  of  P(p)  at  the  double  eigenvalues. 
A  corresponding  peaking  behavior  was  previously  noted   [2,17j 
in  the  experimental  buckling  loads  obtained  by  Kaplan  and 
Fung^'[l].   Although  the  experimental  boundary  conditions 
probably  differ  from  those  of  Problem  B,  this  phenomenon  may 
be  related  to  the  occurrence  of  a  double  eigenvalue. 

Since  the  Jacobian  in  (3.8)  vanishes  if  p  =  P   and 

•^   '  mn 

P  =  P  ,  the  procedure  given  in  Section  3  must  be  modified  to 

investigate  this  case.   The  modifications  are  essentially 

extensions  of  the  methods  used  in  the  Poincare  theory  v/hen 

the  variational  problem  possesses  periodic  solutions,  see 

fl8] .   A  new  parameter  K  is  introduced  and  modified  initial 

value  problems,  ^t   ,  for  the  functions  y   (x:P,e,5,K)  and 

mn  •'  mn  ^  »  '  '  '  ' 

2__  (x;P,  e,cO,K)  are  defined  by 


For  example,  the  experimental  peaks  occur  at  o?-^20,  56  and 
^12  ~^*^'  r'z^^SO.S    .  ' 


•^'  \  •• 


(',)h   'io    er :;;:•.:;;.:;■!    Zy.r.oS.    o-'::.    t5T3i*!T 


£.T 


*:> 


o-i-aOj    '/I       r;jj;ii^,;/..:(.;;'.-    ;},  lb    bx    r-j-    .-  ,    ,      •-■■.    -■    -     JA 


,  f-. 


J.  ■;•;■!,  J,  I 


, -iju  j'evrio  :  .^  o    c-I':/,' oc    ■;■   \-j    tipf!0'T'ix;:ioo    or,:,:    o:t    fjC'i'Ien    ad 

tvfja    ,3noioi:;Xo;i   oiboi-ioa   3c>&ao?;eocj  iueioo'iq   Isfioxrfijl^isv    sric? 

Ijsictini   bBllcbc:.:  bits  feaojjbo'-icfnx   ei    4  'iG.te.u^i^^sq   x^>9n   A      •  i'^il] 

bne    (:U'^  ^3t^'«>0  _^."i   e no. ^Jora;'!   era   •lol    »„,'^   ,eirtorq'o'-;.c3    guIbw' 


'.■JasmxTe 
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:-^v  z   +fcv^™^) 


•'mn    mn''inn   mn  mn     mn   ''inn  mn   '' 


^^"^'  mn    mn  mn     •  mn   "mn 


Y   =0,   y   =1,   z   =0,   z   =6,  forx  =  0 
"'mn    '   ''mn    '    mn    '    mn    ' 


V/e  can  show,  as  in  the  preceding  section  that  the  unique 

solution  of  '^       is  an  analytic  function  of  the  parameters 
mn 

P,  e,  ..^^  and  K.   The  parameter  values  are  to  be  determined 

such  that  the  solutions  of  '^^       satisfy  the  boundary  conditions. 

mn  ''  o  t 


(1^.5)        y^„(liP,e,d,K)    =   z      (l;P,e,g,t)    =   0    . 
inn  liin. 


For  the    special   parameter   values    e   =   >^   =  0   and  P   =  P      =   P 

m    n 

the  solutions  of  -'^  are  y^(x;P^(  P^)  ,0,6^,0)  =  y^^^^xjS^), 

^im(^'^m<^'mn^'°-'L'°)  =  ^^"^^  (^>^°J '  ^^^  ^^^^  ^^^^^^^  (^'^^ 

for  all  bounded  and  real  .^   *   The  implicit  function  theorem 

mn 

is  used  to  determine  other  roots  of  (I4..5)  near  each  root 

rP„,(>-,„v,)  »0,'^„,ol  .   Proceeding  as  in  Section  3  and  using 
■-  m  mn     mn  -10  ->         o 

the  appropriate  variational  problems  we  can  show  that  the 
Jacobian  of  (l+.S)  with  respect  to  P  and  t  and  evaluated  at 

each  root  does  not  vanish  for  all  ^   4  fy     ,  ^    *      Thus  there 

mn  '   'mn  nm 

are  solutions  of  (i^..5)» 


(^.6)     ,.  P  =  P^(e,.^)  .  L=   fc^(e,.s)  , 


61 


,r.-,,\  c     ;„ 


'.■     -•      ,,      .1,0...      ,  f  ;    --      _   _ii  >  ',.'  ,^      .»  J....      •        _      V',  0 '. .' 


it  .T 

^'flit!  >.-  •       »;<;•;■,-••       •  V-;i,f      '■jtj      «  Hii:  '  II;.; 

weioed*   rtoictor.jjl   iioiini'.l    Q.dT     «      c  Ibe-t  ;jnB  bo.h-"joc.   IIjs  -lol 
•?oo'i  riofio   •^•i.:3;;    rc^.ii)    lo   e.'iooT   T.-irljO   tsruitxt&ctvTiL   Oo   fioai..-   ei 


-i  -- 
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which  are  analytic  functions  of  e  and  6  in  some  sufficiently 

small  neighborhood  of  each  root  [p  (/-'  ),0,b   .o]  ,  for  v/hich 

*-  m  '  mn  mn     ^ 

^mn  ^  ^mn'    ^nm  ^^^    satisfy   the    conditions 


(1^.7)        P      (0,6°    )    =  P      =  P      ,      K      (0,5°    )    =  0    . 
'■^  mn'mn  m  n*        iTin^'  nnn 


Hence   the   raodif  ied  problems    (l+.Lj.)  ,    (14-.5)    have    solutions 

7/      \(x:e,6)    ^  Y      (x:P      (e,'^),e,^,  fc      (e,r-)'), 
•'(mn)       *    '  ''mn^    *    mn^    >-'}    f    >   j^^\    t<i)^t 

^(mn)^^'^'^^  -  ^mn^^'^mn^^'^^'^'^^'Sm^^*^^^  v;hich  are  analytic 
in  e  and  6  for  sufficiently  small  lei  and  \S--'^°    I. 

Solutions  of  the  modified  problems  are  solutions  of 
the  original  problem  (3.2),  (3.3)  and  (3.1^)  if  6  and  e  are 
chosen  so  that  the  bifurcation  equations  [l8j 

(i^.8)  ^mn^^*^^  "  °  * 

are    satisfied.     Thus    if    ([(..6)    have    solutions   6=  -'^      (s)    \>ihich 

mn 

satisfy  the  conditions 


then 


f^(x;e)  =  ex-ly(^^(xje,&^(8))  , 


(i^.lO) 

g^(x;e)  =  e-^"^y(^)(x;e,<5^(e))  , 

are  solutions  of  Problem  B  with  f  =  P   • 

'     mn 

To  investigate  the  solutions  of  (I4..8)  we  employ  the 
analyticity  of  y,  2,  P  and  K   and  expand  them  in  the  forms; 


'ci 


•   ;'lr;;^i:>i'ilU3    cr^os    fti  c    bc-fi    '    lo   ^r:-::  il-cn^l    i>s::J[xiaciB    'tis   rioijivj 


c. 


f--:     -  . .  \ 


7;::t',:i.yar    .:.'^/i  ;'oxr:w   ii''^,^^)  ,,.y  ,'-^''  ^^'^  ^'-'^  c^/'^^O  r..f^   ^    '^"^r^  ■ '■•' (niu. 


jji...      I-  ; 


i^'     "[^.:0    -'iJAeioll'li^a   'icih    biU',   3    r;x 


r,;..r    -i.    ^!  ;■      (:!.f   'i     Gil..      >. -J  j  •.  j      j'S.t)     liliiOO'XCj     X :' n/.jj  I'lO     SHo 


HvU"'''^ 


(01,^) 


.    ((3),^^.^^>^)/rm.^^^^'    "    ^=^^-^.1.^ 


nni 


(rm-f) 


.     n  «  ■:\   ricfjw  a  MQldot'i  lo  enp-rd-p- lor.    bib 
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ik'll) 


K — X 

mn    m  -^ — ^  k    ^^ '         '      'mn     r; — =-  k    ^    ' 
K=l  k=l 


where  we  have  used  the  last  of  ([|..7)«   The  series  in  (I4..II) 

converge  in  sorae  interval  lei  <  e   .   The  coefficients  in  the 

°  '  '     mn 

expansions  are  determined  by  substituting  ([|..ll)  into  (Lj-.Li.) 
and  (L|.,5J  •   The  linear  boundary  value  problems  for  y-j^^  and 
z,    ,  obtained  from  the  coefficients  of  e,  possess  solutions, 
if  and  only  if  the  inhomogeneous  terms  satisfy  the  appropriate 

orthogonality  condition.   This  condition  determines  ^;^' 

,  -,(mn) 
and  P-][    as 


(U.12) 


where. 


1    ^^^         ^^mn  -nm^   ^  mn^mn  'nm^nm^  ' 


^1    ^    '         ^  m   n^   ^^mn  ^nm^ 


2  ,s  ■■■ 

Amr.  = -ran  mn  ;j^-2l   a^  +2"^([i+r.2  )(j)   a   a  +(^~2+l)i  ^     ^ 

inn     J2/^^  X  /^  rm^mm  mn     ^^  nm'^mnmnnra  *■  nra  ^  ^nra  nn  '  * 
mn  2^  m' 

(1+.13) 


^mn  =  ■   jf('^m^)Jl("n^)^  ' 
0 


;■■; 


CO 


I~iL 


CO 


fnn-) 


=    '■'■■'  ^^  '■''^h.j:iti)'i 


(X'     -.  -  v\  , 

1  \    --■'/ 


CO 


^i.o>  ■■ 


/■;;■■;•' 


:-;)(-^^): 
V    i/ 


CO 


(Ii.il) 


(Jl.-ji)    i)-;    zerxoz    v^JT      »C^,.ii)    'to   :jo;;I    or.:t    L;^c-.;J    :'\rcri    eu   O'Xgtriv? 


sr;i 


;rl'i:;:>J'3i;  i'..  oC  ;•    oii 


1  .-■  I 


'TO;l;:i     ■:).;K;-t^    ni:     9^;'tOVJ-tOO 


■'.;'.-:]}    c:!--;;    (rr^i'j    ;nn.r.i?M  J  r) -j  ;;:.ih    '/,,'    usrufij'^e^tsb    '^o:;j   ^noienaqxo 


?    ,';r^\ 


^^"'^yi   --!ol    ^ncx-lyici    ^:  I^v  -,;'i:oi:acfOd   •lii-^n.M    oriT      ,  lc:c.;i)    hxi; 


,  BnOlJ-UlCT. 


S.'^.sx'iqoT:.;!; 


ui    ajJ^:crf;--;o:,;oxiriJ;    '■n-^    1x    rluo   bnB   li 


'^"ii    F;':r:J:i;nt' /:.  J:    i..oj.;l  i-rifiOO    e  r:JT      .  ;i'^  cjr.r5noo   vi  .i.l.e<i"ojiOi.':i'to 


{ri::,i) 


r^   b:- 


(ni:;)  , 


(i;x.4) 


t&;?onw 


nui  • 


■■j    nr; 


:.S 


S'.i   'S  m.i 

(tl.-li) 


^     T,  '  ir    '  ' )  ^T, 
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Using  the  last  of  (l+.ll),  the  bifurcation  equations 
(l4.,8)  reduce  to 


(I^.IL,)         ^:/Cj("^^6)e^"^  =  0  . 


If  (U.lii)  have  continuous  solutions  5=6  (e)  with 

mn   ' 

^mn^'^^    ~   "^mn  ^^^^  ^^^   must  satisfy  the  bifurcation  conditions 


(Uci5)  ^I'^^O  =  °  • 


Conversely,  if  S^^^   satisfy  the  bifurcation  conditions  (lj-.l5) 
and  if 


,  (mn) ,  o  . 
(I;. 16)  .-_L-__JS2-  ^  0 


then  the  implicit  function  theorem  is  applicable  to  (ii.lL;) 
and  it  implies  that  the  bifurcation  equations  can  be  solved 

in  some  sufficiently  small  neighborhood  of  3=^°  ,  e  =  0. 

°  -Ton 

The  solutions  ^  =  5^^^^^)    are  analytic  functions  for  suffi- 
ciently small  e  and  satisfy  6     (0)  =  f°    . 

J      mn^      mn 

Inserting  (l|.12)  and  (I4..I3)  into  (1+.15)  and  (I^-.lG)  and 
defining  ^^^  by 

a   (r°  ) 
'  nxa  "Vm^ 

the  bifurcation  conditions  reduce  to  the  cubics: 


^'"'''l-.-^''  {lli..^) 


L-A 


;,■•.;   .-i' -i--,  >(■•",  ;•,-;     .-sri-j-    --■■•,■■  '■■+.;!■.,■     -'-    .-' ,-       ''-'.    r---.     V       '■ '!,    ,-     s' f") ';       ,V 


•"!     ..     i'     ^.-, '  -''■■'   ^  '"J  '■     ,!^ 


;;j.Iu.'0:->   ac.!-:;.--oirjl.i:d    -^r^  y'lt^jy.^    ,.,"      ',11    -\. i^;;;'Tovi:c  0 


'i -/  r     '■' 


(..'l,l')    c:t   o.Cc)irai:I•:;0^^   cii   rfe'ioocict   r!oI.:^::^OJ'^    :M:jJ:j!(iFii    oii:i   cy^rdi 

(          J )        -D 
—    ~  ^  \V-!-*lJy 

Vifn'^  'asi 
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(k.lS)      ^+E   C^+G  E   ^   +G   =0 


and  the  "solvability  conditions"  {l].,l6)    reduce  to: 


/=•''■. 3 ,  r/;-„\3 


C^-l"  4*2(r;,i  ^inHUr)  ^mn-°.mEnm  4-2°mn«»«-(rrj  °™  ^  0  . 


'  n  1  , 


Here  v>fe   have   used   the   notation, 


""        (^-4)^.n^  ^^^  ^  ^  ^2(-J      H 


(^.20) 


p    (1+6^    )    Jo('-    )    'i 

Thus  we  conclude  from  (l4.,l8)  that  for  each  m  and  n  there  is 

at  least  one  real  S   that  satisfies  the  bifurcation  condition 

mn 

and  /  „  =  ^^„,  ^ are  not  roots  of  (I4..I8), 

mn    mn   nm 

If  the  explicit  values  of  the  integrals  k>      >    defined 

•^     ^mn' 

in  (l|»13)  are  known  then  the  roots  of  (li.,l8)  can  be  determined 
using  standard  formulae.   These  integrals  were  numerically 
determined  using  Simpson's  Rule  with  500  mesh  points'"  for  m 
and  n  in  the  range  1  <  m,  n  <  2S*      The  Bessel  functions  were 
evaluated  using  the  procedure  of  [19]  for  arguments  between 
zero  and  eight  and  appropriate  asymptotic  formulae  were  used 


Several  of  the  integrals  v;ere  evaluated  using  1000  mesh 
points.   No  significant  difference  was  observed  with  the  re- 
sults of  the  500  point  mesh.   Double  precision  arithmetic  was 
employed  and  all  calculations  were  preformed  on  the  IBM  70914. 
computer  at  the  AEG  Computing  and  Applied  Mathematics  Center 
of  the  Courant  Institute  of  Mathematical  Sciences.   The  author 
is  indebted  to  Dr.  F.  Bauer  for  conducting  the  computation. 


0  cr      i:~     3    ,2;    c-    ■■'\,  .1+  ,^:^  (6i,j) 


:o:t    oOiJhcT    (oI*4i)    ^'r^nclJ Ibnco  v:>±£i; Jrvloo"   3r!:f    biso 


,rc J- :)•/:; rj-rri    a/;;!"   i:sLi;j    o-/ijr;    oi"   otoH 


;'OX?.ii;)iiO::j   nn.i;ctiiOi;-!l^id    ©rid-   s'ii:'!}-;  i:'jj;3   i.?:;j    ..^ ,/!;   .[r'l-i:   on")   ^tsi^O-i    is 

OT  not   "'»dT!.i:oq  rlaoxil  COH  ii':'''i.H  -dlijH  & 'acr.ci;;:.:;^  3fu;c:i  b9nxjK3T:o.t3b 

STOW  enoLtora-A   I©t,a9G   erJT     *ciS   >  u   ,.ni  >   X    ■?'£-iz?'x   odd  as.  a  bn& 

nsawcted  ad'noxitt;,:>'XG  10I    [^'ij    "iy  eoubaobiq   odJ   -^niaij-  &s;tijtiijav3 

i^jf...  .-v.i    ,    I       I. ..u-.^.-^ ■■■■„,. I,   I,.-,  .M- ; ^i_!.»— i^  "■"— -■' ^ -"T 

•••-'(    0001    ^nizu    f-       :---•• -      ■     .;-    -.     --...!-     .  .;-}    -i^    iB'I9V9?, 

'    f^c^iw  £)ov^«p>;-  . 'i   oil     iSifrii     i 

Cilia'   ^3    • 
t  ■  ■  .   jiiioo  Qdi  2«i;js,...;  ,  .....,■•.    -  ,    ,    -I  oi  j^,  .• .  .  .,..^.    .^ 
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for  larger  arguments,  lie   found  that  all  the  roots  of  ([|.,l8), 
for  m  and  n  in  the  rangei"  1  <  m,  n  <  25>  satisfied  (L1..I9) 
well  within  the  accuracy  of  the  coraputation.   In  fact,  the 
left  side  of  (I4.0I9)  was  usually  quite  large.   Hence  we  have 
shown  that  there  is  at  loast  one  solution  of  Problem  B  for  P 
in  some  sufficiently  small  neighborhood  of  each  double 
eigenvalue  with  m  and  n  in  the  above  range.  We  conjecture 
that  there  are  solutions  near  each  double  eigenvalue  for  all 
m,n=l,2,...  .   V/e  also  find  that  for  the  following  indiciea 
(I].. 18)  has  three  real  roots: 

1)  n=m+l,  3<n<25,  n  odd, 

2)  n=m+3,  ll+<n<2St    n  odd, 

3)  n=2m+l,  2<n<25, 
L|.)  n=2m+3,  6<n<25, 
5)  n=2m+5,  ^l"^!^^' 

Thus  in  a  sufficiently  small  neighborhood  of  the  double 
eigenvalues  with  the  above  indicies  there  are  three  solutions 

of  Problem  B.  Hence,  as  O — ^  ^        the  m-th  and  n-th  branches 

I  mn 

coalesce  and  when  P-  p        one  solution  may  be  "destroyed"  or 

1    ■  mn 

a  third  solution  may  be  "created"  depending  on  the  values  of 
the  indicies  m  and  n. 


'  By  symmetry  we  need  only  consider  n  >  m. 


iV£i':l    sw   '^omU      /5ri-.:Bi    :rji:/.;'    ^;;II;^^•;sl;;   a;^'-    {Vi^«-4)    'lo   el^cc    ^l^i: 


T    'X  .■!    &    ;r.j,[dO' 


no  .i:  ■.!■;.!  Cari    t^no   yC/:..f    o'i;.   t,  I    oirtii;'    ji;fid'  riwoiiE 


n^xorib'-c:    ;*:■ ;  viol.' o?:    e;;.]-   i;o;;    c.^;r;i•  t/:-;:';;   o-;i^    -j"i      ,    .  » .  ,StX-n,:^i 


bi'O   :-i    ,.i?'-- 


,;:-h"x5     (j; 


^^^ 


t  -       'i  :.>— ..  1 


C  ^r.S-r: 


^  ncrr 


'to   B9JJl.cy"   erf  J   //.o  an."  j    ''beciBS'-.rrj"    so   ^J"^-'  J^^o 


B 


-■-bnx   3ri;t 


2k 


5,  The  Intermediate  Bucklinp;  Load 

It  is  convenient  to  reformulate  Problem  B  to  establish 
the  existence  of  the  intermediate  buckling  load  P..,  or 
equivalently  Xj^  =  PPm»  ^^^d  to  obtain  upper  and  lower  bounds 
on  its  magnitude.   Equation  (2.1b)  is  integrated  and  the 
second  of  (2.3)  is  used  to  obtain 

X 

(5.1)  g'(x)  =  -f:x-^/   rf2(x^)  +2f(x^)]x^dx^  . 

-0 

Then,  using  (2.l4.b),  g(x)  is  given  as  a  functional  of  f(x)  by 

X 

(5.2)  g(x)  =  ;  g'(x, )dx   . 

-'1 

The  boundary  value  problem,  B,  is  reformulated  as  Problem  B 
as  follows:   To  find  a  function  f(x)  which  possesses  a  contin- 
uous second  derivative  and  satisfies  the  differential  equation 
(2.1a)  and  the  bo\indary  conditions 

(5.3)  f'(0)  =  f(l)  =  0  . 

The  function  g(x)  in  (2.1a)  is  defined  by  (5.1)  and  iS»2) , 
These  equations  imply  that  g(x)  satisfies  (2,1b)  and 
g  (0)  =  g(l)  =  0.   This  statement  of  the  boundary  value 
problem  is  similar  to  the  one  introduced  in  [20]  in  a  study 
of  the  buckling  of  circular  plates. 

The  "energy"  functional  V  is  defined  as, 

1 

iS.k)    V[f(x)>X,f]  ^  ;'"[2(f'^(x)-Xf2(x))  +g'^(x)]x^dx  , 

-0 


.0^ 


•; X.  v..i-. '  ^BiM^i'SL 2jAlk!^-:"iS iixL^.ilX  •  ^ 


,     A./     ,.^ 


k '"- ; 


vd    (x;l    "TO   ifino.i. J Dfi:/'' 


«•-''.' 


n  [  ,.  X ;    ,s 


( i )  ■! 


-   (^0 


U'.c) 


i  Ji: 


ii^    {.r«'i^)    iu   bioiii":.i^k 


Pil'^j   frr:/c!3i-«/i:    -ynvl 


,    ;t&^:MUr'3'^  (^^l^'i^  -(xJ-^'ajSj     ^3.    !  :,A4l;i:)^3V    (m,?) 
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where  g  (x)  is  considered  as  a  functional  of  f(x)  defined  by 
(5.1).   Thus  V  is  a  functional  of  f(x)  only  and  it  is  pro- 
portional to  the  difference  between  the  potential  energies  of 
a  buckled  and  the  unbuckled  state. 

The  relationship  bet'.ireen  B   and  V  is  obtained  by  first 
defining  a  class  of  functions.   f(x)  is  contained  in  A  or  is 
an  A-function  if  in  the  interval  0  <  x  <  1  it  is  continuous 
and  satisfies  (5.3) »  f  (x)  is  an  Lp  function  and  all  integrals 
in  (5.1)  »  {5*2)    and  {Sci\.]    exists   Then  it  is  easy  to  show 
that  if  r(x)  is  a  solution  of  b'  it  makes  V  stationary  with 
respect  to  all  A-functions,  The  converse  of  this  result  can 
be  proved  using  the  raethods  outlined  in  [ZOj  •   The  converse 
states  that  if  f(x)  makes  V  stationary  (or  minimizes  V)  then 
f"(x)  is  continuous  and  f(x)  solves  B.   Here  f  (x)(r.  A  is  aaid 
to  minimize  V  for  a  fixed  X  and  p  if  Vff  (x)]  <  Vff (x)!  for 
all  f  (x)  -^  A. 

The  result  concerning  the  existence  of  X,,  is  nov^r  stated 
as  the 

THEOREII,   If  for  every  finite  P  and  X  there  is  an 
A-function  which  minimizes  V,  then  \iir)    exists  and  it  is  in 

the  interval  '^  <  >^j,j(  r )  <  2i  -  ^  Z^  ^  ^  >   vJhere  P  (p )  is  defined  in 
(^..3)  ,  and  ^'  is  the  first  zero  of  J^  (x)  =  0. 

The  theorem  is  a  direct  consequence  of  the  lemmas  given 
below  and  the  following  inequality,  which  is  easily  demon- 
strated using  classical  methods  [^2lj  , 


zi^'Ti  • '  '^'^'■s    l.i«   b.ii^   ::'-!^"^.  .;' 


V:"'  '■'  p.     :i 


;">    :;  s. 


>.'  .1. 


r-  n 


--    '..lib    :..-:^'     .;J    i.ii'i:^-::-j-:fO'; 
.^    -/:•   :-..:■    li    iu:.•::^.^:.■i..';;,~A   r^: 
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.2^2 


iS.S)    I    ff'(x)  -oJf'^(x)Jx^dx  >  0  ,   for  all  f(x)^  A  . 

p 
Here  ^-',  is  also  the  lovrest  eigenvalue  of  the  linearized 

buckling  problem  for  the  sy/imetric  deformations  of  radially 

compressed  and  clamped  circular  plates: 

(5o6)    Gf(x)  +Af(x)  =  0  ,   f'(0)  =  f(l)  =  0  . 

V/e  call  (5»6)  the  "equivalent"  circular  plate  problem.   The 
inequality  (^oS)  and  the  functional  V  given  in  (5»h-)  and  (5.1) 
immediately  yield 

LEI-MA  1.   If  >^  <  -->?  then  only  f  (x)  =  0  minimizes  V, 
The  form  of  the  functional  V  and  the  properties  of  the 

eigenvalues  (3«5a)  and  eigenvectors  (3c7)  of  the  linear  shell 

buckling  problem  yield 


LEIIIIA  2.   If  X  =  X(p)  then  there  exist  for  each  p  >  0 
a  function  f(x)£A  such  that  V[f(x)J  <  0. 

The  proof  of  this  lemma  is  obtained  by  a  simple  calcula- 
tion of  vff  (x)l  where  f  (x)  =  B  J,  (<^-'  x)/x  are  eigenfunctions 
'-  n   -■        n       n  1   n  ' 

of  the  linear  theory  and  B   are  constants  restricted  to  lie 

in  the  ranges  0  >  B  >  -2B  /a  ,  n  =  1,2,.,.  .   The  constants 

n      nn 

a  and  6^  are  defined  by 
n      n  " 


a2  =  .2  / 
n  ^  ' 


Pn  ^  2p2  / 
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Lemma  1  implies  that  if  A  =  X  then  ^.l.b.  V[f(x)]  <  0.   The 
remaining  two  lemmas  show  that  the  set  of  X ' s  is  divided  into 
tv>fo  disjoint  sets:   those  X's  for  which  the  minimum  of  V  is 
zero  and  those  for  which  it  is  negative, 

LEM>'iA  3-   If  ^(^)  =  0  minimizes  V  for  X  =  X  ,  then  only 
f (x)  =0  minimizes  V  for  all  X  <  X  , 

Proof:!*  Ue  deduce  directly  from  the  form  of  V  that  if 
^  <  ^-^>  V  ^  0  for  all  A-functions.   Suppose  there  is  a 
X  =  X"""  <  X   and  an  A-f unction  f'"(x)  ^  0  which  minimizes  V, 
i.e.  V[f'"'(x)  ;X'"]  =  0.   Then  if  X  is  in  the  range  X'"'  <  X  <  X  , 
V[r"iXj  <  Vff"'''>X'''J  =  0.   This  is  in  contradiction  to  the  non- 
negative  property  of  V  for  X  <  X   and  the  proof  of  the  lemma 
is  complete, 

LEI-IMA  I4..   If  the  minimum  of  V  is  nep;ative  for  X  =  X 
then  it  is  negative  for  all  X  >  X  . 

ii _      Q 

Proof:   By  contradiction  using  Lemma  3» 
The  theorem  is  then  proved  if  we  define  ^T,r(  "')  as  the 
l.u.b.  of  those  X  for  which  f (x)  =  0  minimizes  V.   Thus  for 

each  f   there  is  a  X  =  X„  such  that  for  all  X  >  X.,  there  are 

*  ri  M 

buckled  states  with  less  energy  than  the  unbuckled  state  and 
for  X  <  X.^,  all  buckled  states  have  greater  energy  than  the 
unbuckled  state. 

The  intermediate  load  can  also  be  characterized  by  the 

'  This  form  of  the  proof  v/as  suggested  by  Dr.  II.  Nei-;man. 
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minimuin  property. 


(5.7)    \A.)    =  i^in 

^^  '     f  (x)  ?^  0  €A 


■'  (2f '2^  g'2)x3dx 


I 

"o 


J 


where  g' (x)  is  defined  in  (5.1).   Upper  bounds  for  X^^  are 
obtained  by  selecting  trial  A-functions  to  make  the  quotient 
in  (5.7)  as  small  as  possible.   For  example,  functions  of  the 
form  f(x)  =  pF(x)  are  considered,  where  P(x)  -A  is  a  specified 
fimction  and  (3  is  an  arbitrary  constant.   The  quotient  is 
then  a  function  of  p  only  and  {3  is  determined  to  minimize  it. 

Some  of  the  trial  functions  that  were  used  and  the  resulting 

—  32 
upper  bounds,  X,,  .(  )  =  a.  +b.lO   r^  ,  are  shown'""  in  Table  I, 

1 1  j  X  X     X 


Table  I 


i 

F(x) 

^M,i  =  ^i 

+  b.io-V 

^ 

^i 

1 

1-x 

15.0 

6.5166 

2 

l-x2 

16.0 

3.7037 

3 

l-x3 

17.50 

2.!+l+5l 

U 

l-x^ 

19.20 

1,7166 

5 

l+x-2x2 

16.7576 

2.7636 

6 

2(2-x-x2) 

15.5051^. 

t^.6055 

The  calculations  i;ere  performed  by  M,  Szeto. 
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The  results'"  are  suranarized  in  Pig,  3  v;here  v;e  also  show 
accurate  approximations  of  A,,  obtained  from  a  numerical 
solution  of  Problem  B.   Details  of  these  and  other  calcula- 
tions will  be  reported  elsewhere  [22]  . 


6.  Other  Bifurcation  Problems 

Analogous  results  can  be  obtained  for  other  bifurcation 
buckling  problems.   For  example,  for  the  boundary  value 
problem  consisting  of  (2.1),  (2,3)  and  (2e5)  which  we  call 
Problem  3-,  ,  we  can  show  that 

vjhere  X  is  now  the  minimura  eigenvalue  of  the  linearized  shell 
buckling  theory  using  the  boundary  conditions  (2.5)  in  place 
of  {2,k)  '      Upper  bounds  for  ^^^(p)  are  obtained  from  a  formula 
similar  to  (5.7)  .   These  results  are  ^•^aphed  in  Pig.  5  'vith 
the  predictions  of  ?v„  obtained  from  a  numerical  solv-.tlon 


of  the  boundary  value  problem   [22] 
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Captions  for  Fipiures 

Figure  1.   Sketch  of  conjectured  load  deformation  curve  for 
bifurcation  buckling.   Here  D  is  a  representative 
deformation,  e.g.  the  normal  displacement  of  the 
cap's  center. 

Figure  2.   Shell  geometry. 

Figure  3«   Intermediate  buckling  loads  for  Problems  B  and  B'"'. 

The  upper  bound  curves  are  essentially  the  envelopes 
of  the  curves  given  in  Tables  I  and  II. 
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Attn:  Op  91  (1) 

Op  03EG  (1) 


Commandant,  Marine  Corps 
Headquarters,  U3MC 
Washington  25,  D.C, 


(1) 


Office  of  Chief  Signal  Officer 
Dept,  of  Army,  Washington  25,  D.C, 
Attn:  Engin.  and  Techn.  Div.     (1) 
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Chief,  Buream  of  Ships 

Dept.  of  Navy 

Washington  25,  D.C. 

Attn:  Code  106 

(1) 

Code  312 

(5) 

Code  320 

(1) 

Code  370 

(1) 

Code  375 

(1) 

Code  U20 

(1) 

Code  U21 

(1) 

Code  U23 

(2) 

Code  U2$ 

(1) 

Code  UiiO 

(1) 

Code  UIi2 

(2) 

Code  Uii3 

(1) 

Code  525 

(1) 

Code  633 

(1) 

Cliief,  Buream  of  Aeronautics 

Dept.  of  Navy 

VJashington  25,  D.O. 

Attn:  AE-U 

(1) 

AV-3U 

(1) 

AD 

(1) 

AD-2 

(1) 

RS-7 

(1) 

RS-8 

(1) 

SI 

(1) 

AER-126 

(1) 

Chief,  .3ureani  of  Ordnance 

Dept.  of  l--avy 

Washington  25,  D.C. 

Attn:  Ad3 

(1) 

Re 

(1) 

ReS 

(1) 

ReU 

(1) 

ReS5 

(1) 

ReSl 

(1) 

Ren 

(1) 

Spec.  Proj.  Office,  3ur.  Ord. 

Dept.  of  Navy 

V/ashington  25,  D.C. 

Attn:  iiissile  ^r. 

(2) 

Chief,  Bur.  Yards  and  Docks 

Jept.  of  Navy 

Washington  25,  D.C. 

Attn:  Code  D-202  (1) 

Code  D-202. 3  (1) 

Code  220  (1) 

Code  D-222  (1) 

Code  D-IilOC  (1) 

Code  D-lihO  (1) 

Code  D-500  (1) 


Commanding  Officer  and  Director 

David  Taylor  Model  3asin 

Washington  7,  D.C. 

Attn:  Code  liiO  (1) 

Code  600  (1) 

Code  700  (1) 

Code  720  (1) 

Code  725  (1) 

Code  731  (1) 

Code  7U0  (2) 

CO,  U.S.  Naval  Ordnance  Lah. 
White  Oak,  Maryland 

Attn:  Techn.  Library  (2) 

Techn.  Eval  Dept,         (l) 

Director,  materials  Lab. 

N.Y.  Naval  Shipyard 

Brooklyn  1,  N.Y.  (D 

CO,  Portsmouth  Naval  Shipyard 
Portsmouth,  New  Hampshire        (2) 

CO,  uare  Island  Nav.  Shipyard 
Vallejo,  California  (2) 

CO  and  Director 

U.S.  Nav,  '■electron.  Lab. 

San  Diego  52,  California         (1) 

Officer-in-Charge 
Nav.  Civ.  Engine  Res. 
and  Eval.  Lab. 
U.S.  Nav.  Jonstr,  Battal.  Center 
Port  Hueneme,  California         (2) 

Dir.,  Nav.  Air  Experimental  Sta, 
Nav,  Air  Mat.  Center,  Nav.  Base 
Philadelphia  12,  Penna^ 
Attn:  Materials  Lab.  (1) 

Structures  Lab.  (l) 

Officer-in-Charge 

Underwater  Lxplos.  Res.  Div, 

Norfolk  waval  Shipyard 

Portsmouth,  Virginia 

Attn:  Dr.  A.H.  Keil  (2) 

CO,  U.S.  Nav.  Proving  Ground 
Dahlgren.  Virginia  (1) 

Supr.  of  Shipbuilding 

USN  and  Nav.  Inspec.  of  Ordnance 

General  DjT.amics  Corp., 

Electr.  Boat  Ldv< 

Groton,  Connecticut  (1) 
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Supr.  of  Shipbuilding 

USN  and  Nav,  Inspec,  of  Ordnance 

Newport  News  Shipbuilding 

and  Dry  Dock  Co, 
Newport  News,  Virginia         (1) 

Supr,  of  Shipbuilding 
USH  and  Nav,  Inspec.  of  Ordnance 
Ingalls  Shipbuilding  Corp< 
Pascagoula,  Mississippi        (1) 


CO,   U.S,  ilav.   Admin.   Unit 
KIT,   Cambridge  39,   iiass. 

Officer-in-Charge 
Postgrad,  School  for 

uaval  Officers 
Webb  Inst,  of  Nav,  Arch, 
Crescent  Jeach  Rd. 
Glen  Cove,  L.I.,  il.Y. 

Supt.,  Nav.  Gun  Factory 
VJashington  25,  D,C. 

Comm.,  Nav,  Ordnance  Test  Sta, 
China  Lake,  California 
Attn:  Physics  Div. 

Mechanics  Div, 

CO,  Nav.  Ordnance  Test  Sta. 
Underwater  Ordnance  Div, 
3202  E.  Foothill  Jlvd, 
Pasadena  8,  California 
Attn:  Struc.  Div, 


(1) 


(1) 


(1) 


(1) 

(1) 


(1) 


CO  and  Director 

U.S,  Nav.  Engin.  txp, 

Annapolis,  Maryland 


station 


(1) 


Supt.  U.S.  Uav,  Postgrad,  School 
Monterey,  California 

Comm.  Marine  Corps  Schools 

Quantico,  Virginia 

Attn:  Dir.,  MC  Dev,  Center 

Comm  Gen.,  USAF 
Washington  25,  D.C. 
Attn:  Res.  and  Dev,  Div, 

CO,  Air  Material  Comriiand 
Wright-Patterson  AF'o,   Ohio 
Attn:     MCHEX-D 

Structures  Div, 

CO,  USAF  Inst,  of  Technology 
VJright-Patterson  kFB,   Ohio 
Attn;  Chief,  Appl,  Mech. 
Group 


(1) 


(1) 


(1) 


(1) 
(1) 


(1) 


Director  of  Intelligence 
Headquarters,  USAF 
VJashington  25,  D.C, 
Attn:   PV  jr.  (Air  Targ,  Div) 

CO,  A.F,  Office  3ci,  Research 
'; Washington  25,  D.C, 
Attn:  Mechanics  Div, 

U,S,  Atomic  Energy  Commission 
Uashington  25,  D,C, 
Attn:  Dir,  of  Research 

Dir,,  Nato  Bur.  of  Standards 
iJas'aington  25,  D.C, 
Attn:   Div,  of  Mechaiiics 
Engin,  Mech,  Sect, 
Aircraft  Structures 

Comm.,  U.S.  Coast  Guard 

1300  £.  St.,  iJV; 

..ashington  25,  D.C, 

Attn:   Chief,  Test  and  Dev,  Div, 


(1) 


(1) 


(2) 


(1) 
(1) 
(1) 


(1) 


U,S.  Maritime  Administration 

General  Admin,  Office  Jldg, 

Lial  G  St.,  NV; 

t/ashington  25,   D.C. 

Attn:  Chief,  Div.  Prelim. 

Design  (1) 

Nat,  Aero,  and  Space  Admin. 

Langley  iissearch  Jenter 

Langley  Field,  Virginia 

Attn:   Structures  Div.  (2) 

Nat,  Aero,  and  Space  Admin, 

1512  H  St.,  Wii 

VJashington  25,   D.C, 

Attn:   Loads  and  Struc.  Div,       (2) 

Director,  Forest  Prod,  Lab. 

Madison,  Wisconsin  (1) 

Federal  Aviation  Agency 

Dept,  of  Commerce 

Wasliington  25,  D,C, 

Attn:   Chief,  Air  Lngin,  Div,     (1) 

Chief,  Air  and 

Equip,  Div,  (l) 

National  Science  Foundation 

1520  H,  St.,  mi 

VJashington,  D.C,  (l) 

General  Dynamic  Corp, 

i^lectr,  Joat  jiv, 

Groton,  Connecticut  (l) 
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National  Academy  of  Sciences 

2101  Constitution  Ave., 

Washington  25,  D.C, 

Attn:  Dir-,  Conn,  on  Ships 

Struc.  Jesign  (1) 

Exec.  Secy.,  Oomn,  on 
Undersea  ..arfare        (1) 

Newport  Hev/s  Shipbuilding  and 

Dry  Dock  Co. 
Newport  News,  Virginia         (1) 

Ingalls  Shipbuilding  Corp, 
Pascagoula,  Mississippi        (1) 

Prof,  L^Tin  S,   3eedle 

Fritz  Engineering  Lab. 

Lehigh  University 

Bethlehem,  Penna.  (1) 

Prof.  R.L.  Bisplinghoff 
Lept.  of  Aero.  Lng^neering 
Massachusetts  Inst,  of  Techn, 
Cambridge  39,  iiassachusetts     (1) 

Prof.  H.  !.  Sleich 

Dept,  of  Civ.  Engineering 

Columbia  University 

New  York  27,  Hew  York  (1) 

Prof.  3. A,  3oley 

Dept  of  Civ.  Engineering 

Columbia  University 

New  York  27,  New  York  (l) 

Dr.  John  F.  Brahtz 

Southern  California  Labs 

Stanford  Research  Institute 

820  Mission  St, 

South  Pasadena,  California      (l) 

Dr.  D.O.  Brush 

Struc.  Dept.  53-13 

Lockheed  Aircraft  Dorp. 

Missile  Syst.  Div. 

Sunnjn/^ale,  California  (1) 

Prof.  B,  Budiansky 

Dept.  of  Hech,  Engineering 

School  Appl.  Sciences 

Harvard  University 

Cambridge  3&,  Massachusetts     (1) 

Prof.  Herbert  Deresiewicz 

Dept,  of  Civ.  Engineering 

Columbia  University 

632  Vi.  125th  St. 

New  York  27,  N.Y.  (1) 


Prof,  D.C,  Drucker,  Chmn, 

Div,  of  Engineering 

Brown  University 

Providence  12,  Rhode  Island      (1) 

Prof.  John  Duberg 

Dept.  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof,  J.  Ericksen 

Mech,  ..ngineering  Deot. 

Jolins  Hopkins  University 

Baltimore  16,  Maryland  (1) 

Prof,  A.J,  Eringen 

Dept.  of  Aero  Engineering 

Purdue  University 

Lafayette,  Indiana  (1) 

Prof,  ir,  Flugge 

Dept.  of  Hech,  Engineering 

Stanford,  California  (1) 

Mr.  M.  Goland,  VP  and  Dir. 

Southv;est  Research  Institute 

8500  Culebra  lid. 

San  Antonio  6,  Texas  (l) 

Prof.  J.N.  Goodier 

Dept  of  hech.  Engineering 

Stanford  University 

Stanford,  California  (l) 

Prof.  L.'!,  Goodman 
i.ngineering  Ixperimental  Sta, 
University  of  Minnesota 
Minneapolis,  Minnesota  (1) 

Prof,  H,  Hetenyi 

The  Technical  Institute 

Northwestern  University 

Evanston,  Illinois  (l) 

Prof.   ?.G.  Hodge 

Dept.   of  Mechaiiics 

Illinois  i^.nst.  of  Technology 

Chicago  16,  Illinois  (1) 

Prof,  N.J.  Hoff,  Head 

Div.  Aeronautical  Engineering 

Stanford  University 

Stanford,  '^alifomia  (1) 

Prof,  W.H,  Hoppmann,  II 

Dept,  of  Mechanics 

Rensselaer  Poljitechnic  Inst, 

Ti'oy,  i>iew  York  (1) 
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Prof.  3ruce  G,  Johnston 

University  Of  Michigan 

Ann  Arbor,  Michigan  (1) 

Prof.  J.  Kempner 
Dept,  of  Aero,  Engineering 
and  Appl,  Mechanics 
Polytechnic  Inst,  of  Brooklyn 
333  Jay  3t. 
Brooklyn  1,   ll.Y.  (1) 

Prof.  H.L.  Langhaar 
Dept.  of  Theoretical  and 
Applied  Mechanics 
University  of  Illinois 
Urbana,  Illinois  (1) 

Prof,  B.J.  Lazan,  Director 
Engineering  .jqjeriiaental  Cts. 
University  of  Minnesota 
Minneapolis  lU,  Minnesota      (1) 

Prof.  L.H.  Lee 

Div.  of  Appl,  Mathematics 

Brown  University 

Providence  12,  Rhode  Island     (l) 

Prof.  George  H.  Lee,  Dir.  of  Ues. 

Rensselaer  Polytechnic  Inst. 

Troy,  N.Y.  (1) 

Mr.  S.  Levy 

GE  Electr.  Research  Lab. 

3198  Chestnut  St. 

Philadelphia,  Penna.  (1) 

Prof.  Paul  Lieber 

Geology  Departnent 

University  of  California 

Berkeley  h,   California         (1) 

Prof.  Joseph  ilarin,  Head. 
Dept.  engineering  Mechanics 
College  of  Lngin.  and  Arch, 
Pennsylvania  State  University 
University  Park,  Penna.        (l) 

Prof.  K.D.  Mindlin 

Dept.  of  Civ.  Engineering 

Columbia  University 

632  W.  125th  3t. 

New  York  27,   Wcw  York 

Prof.  Paul  H,  iJaghdi 

Building  T-7 

College  of  Engineering 

University  of  California 

3erkeley  U,   California  (1) 


5  - 

Prof.  William  A.  Hash 

Dept  of  Engineering  Mechanics 

University  of  Florida 

Gainesville,  Florida  (1) 

Prof.  H.M.  Hewmark,  Head 

I'ept  of  Civ.  Engineering 

University  of  Illinois 

Urbana,  Illinois  (1) 

Prof,  E.  Orowan 
Dept.  of  Mech,  Engineering 
Massachusetts  Institute  of  Techn. 
Cambridge  39,  Massachusetts       (1) 

Prof,  -^ris  Phillips 

Dept,  of  Civ.  Engineering 

15  Prospect  St. 

Yale  University 

New  Haven,  Connecticut.  (1) 

Prof.  vJ.  Prager,  Chmn, 

Phys.  Sci.  Council 

:Jro'.>m  Univc-rsity 

Providence  12,  Rhode  Island       (1) 

Prof.  J.R.M.  itadok 

Dept.  of  Aero  ont^rieering 

and  Appl.  Mechanics 
Polytechnic   Inst,   of  Jrookl;;>Ti 
333  Jay  St. 
Brooklyn  1,   W.Y.  (1) 

Prof.   -.L.  Reiss 

Inst,   of  Piathenatical  Sciences 

New  York  University 

li  V/ashington  Place 

New  York  3,   i!.Y.  (1) 

Prof,  ^.  Reissner 
Dept,  of  Mathematics 
Massachusetts  Inst,  of  Technology 
Cajabridge  39,  Massachusetts       (1) 


Prof.  M.A.  Sadov;3ky 
Dept.  of  Mechanics 
Rennselaer  Polyteclinic  Inst. 
Troy,  New  York 

Prof.  J.vV.  Shaffer 
De;)t.  of  Mech,  Engineering 
Hew  York  University 
University  Heights 
New  York  53,  N.Y. 

Prof.  J.  Stallmeyer 
Dept.  of  Civ.  Engineering 
University  of  Illinois 
Urban,  Illinois 
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Prof.  Eli  Sternberg 
Dept.  of  Mechanics 
Brown  University 
Providence  12,  Rhode  Island 

Prof.  T.Y.  Thomas 
Grad.  Inst.  Math  and  Mech, 
Indiana  University 
Bloomington,  Indiana 

Prof,  3.:^,  Tinoshenko 
School  of  engineering 
Stanford  University 
Stanford,  California 

Prof.  A.S,  Velestos 
Dept  of  Cxv,   Engineering 
University  of  Illinois 
Urbana,  Illinois 

Dr.  W,  v;enk 

Southwest  lie  search  Institute 

8500  Culbera  Rd. 

San  Antonio,  Texas 

Prof,  Dana  Young 

Yale  University 

New  Haven,  Connecticut 

Prof.  R.A.  Di  Taranto 
Dept.  of  1-iech.  .Jngineering 
Drexel  Institute 
32nd  and  Chestnut  Streets 
Philadelphia,  ?enna. 

Mr.  H.K.  Koopman,  Secy, 
Welding  Res,  Council 

EngineerinK  Foundation 
29  '-/.  39th  St, 
New  York  18,  M.Y. 

Prof,  '..alter  R.  Daniels 
School  of  hn^ln.   and  Arciiit, 
Howard  University 
Washington  1,  D.C, 

Comin,,  (Code  753) 
U.S.  Naval  Ordnance  Test  Sta, 
China  Lake  California 
Attn:  Techn,  Library 

Prof.  J.  i  .  Cermak 
Dept,  of  Civ,  Ln(3ineering 
Colorado  State  University 
Fort  Collins,  Colorado 


(1) 


(1) 


(1) 


(1) 


(1) 


(1) 


(1) 


(2) 


(1) 


(1) 


(1) 


Prof,  '..',J,  liall 
Dept,  of  Civ,  Engineering 
University  of  Illinois 
Urban,  Illinois 

Dr.  HjTian  Serb  in 
Design  Integration  Eept. 
Hughes  Aircraft  Co, 
Culver  City,  California 

CommaJider,  V.'ADD 
Wright-Patterson  .-iFB,  Ohio 
Attn:   W;iJRC 

VMWiDS 

Comnanding  Officer 

USNUOEU 

Kirtland  Air  Force  3ase 

Albuquerque, 

Httn: 


Nev7  Mexico 
Code  20 
(Dr.  J.N,  Brennan) 


Legislative  ;-:eference  Sei'^ace 
Library  of  Congress 
Washington  25,  D.C, 
Attn:   Ur.  E.  -enk 


Dr,  A,  Ross 

Aircraft  Nuclear  Propulsion  Lept. 
General  ..lectric  Co, 
Cincinnati  15,  Ohio 

Dr,  F,  Lane 

General  Applied  Science  Labs 
Stewart  and  Merrick  Avenues 
l.estbury,  L.I.,  M.Y. 

Commander  Ldward  Leonard 
Asst,  ilavy  Representative 
KIT  Lincoln  Lab, 
Lexington  73,  Massachusetts 

Prof,  W,  Pohle 
Grad,  Math,  Dept. 
Adelphi  City,  L.I.,  N.Y. 

Dr.  Martin  Goldgerg 
tiesearch  Dept, 
Grunann  Aircraft 
Jethpage,  L.I.,  iI.Y. 

Prof,  3,D,  Larmarajein 
Aero  ^pace  De-t, 
San  Diego  State  University 
San  Diego  15,  California 


(1) 


(1) 


(1) 
(1) 
(1) 
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(1) 


(1) 


(1) 


(1) 


(1) 
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